function raux1 = liouv(nd,rho,omsc,gamsc,cset,coef)
% Action of the scaled Liouville superoperator over rho.
%                  L rho = raux
% WARNING: This subroutine is model and representation dependent!!!

         raux1 = 0.;
          for ii = 1:nd*nd
         if (abs(coef(ii))>0.) 
          auxop = cset(:,:,ii);

      raux1 = raux1 + coef(ii)*auxop'*rho*auxop;
      raux1 = raux1 - coef(ii)*0.5d0*auxop*auxop'*rho;
      raux1 = raux1 - coef(ii)*0.5d0*rho*auxop*auxop';
          end 
          end 
          
       end

